Method and device for controlling angular speed of an electromechanical chain with low damping

ABSTRACT

The invention concerns a method and a device for controlling the angular speed of an electromechanical chain with low damping. Speed control of a load ( 10 ) driven in rotation by an electromechanical chain with low damping, possibly having transmission nonlinearities such as clearances, without measuring said speed, is provided by jointly using an observer ( 18 ) and a robust corrector ( 20 ) of non integer order, of the CRONE type, whereof the physical execution is provided by a corrector of integer order, obtained by approximation in the frequency domain. The device thus provides proper damping of the system in the presence of resonant frequencies and excellent robustness, while enabling a relatively simple adjustment.

[0001] The invention concerns a method of controlling the angular mechanical speed of a load driven in rotation by an electromechanical drive train with little damping, and possibly having non-linearities such as transmission play, in the absence of a measurement of angular speed of the load.

[0002] The invention also concerns a device for implementing this speed control method.

[0003] Throughout the text, the term “load” designates indiscriminately any mechanism that can be driven in rotation by means of an electromagnetic drive train with little damping. Thus, without being limited thereto, the load can be composed of a roller of a rolling mill or paper manufacturing machine, a ship propeller, a railroad drive axle, etc.

[0004] Moreover, the expression “electromagnetic drive train” designates an assembly including a torque actuator formed by a power converter and an electric motor of any type, as well as a drive train through which this actuator drives the load.

[0005] Finally, the expression “with little damping” means that the electromechanical drive train has mechanical resonance frequencies. These frequencies are defined by the drive train's parameters, such as inertia and rigidity. In the context of the invention, these parameters are not well known or vary during operation, which results in uncertainty or variation of the resonance frequencies of the electromagnetic drive train and of the associated amplitudes.

[0006] Diagramatically represented in FIG. 1 of the attached drawings is a known type of closed loop servo system providing the control of the speed of rotation of a load 1, driven by a torque actuator 2 through a drive train 3. The torque actuator 2 includes an electric motor 4 coupled to a power converter 5.

[0007] The control of the speed of the motor 4 is provided by a controller 6 that receives a reference signal Ω_(m)* representative of the desired speed of the motor 4. When there is no measurement of the speed of the load 1, the controller 6 also receives a signal Ω_(m) representative of the actual speed of the motor 4. In response, the controller 6 issues a torque reference signal Γ*, which controls the torque actuator 2.

[0008] Among existing controllers, for example, are PI (Proportional Integral), PID (Proportional Integral Derivative), H_(∞), QFT (Quantitative Feedback Theory), and RST.

[0009] In the frequency domain, and using Laplace transforms, it is known that the response of a controller is characterized by a rational transfer function C(s), which corresponds to the ratio of two polynomials N(s)/D(s), where s designates the complex operational variable of the system in question. The transfer function of the open loop system that includes the controller is represented in the Black, Nyquist, or Bode plan, by a curve called “frequential pattern,” which shows the changes in gain and phase difference between the system's input and output, as a function of the real part ω of the operational variable s.

[0010] The different controllers are distinguished from each other particularly by the degree of each of the polynomials N(s) and D(s), that is, by their number of zeros (roots of the numerator) and poles (roots of the denominator). When the system is high order, with little damping, that is, when it has a high number of poles and zeros, only a high order controller makes it possible to meet the objectives of robustness in performance and stability. However, this also results in increasing the number of coefficients of the controller's polynomials, and this consequently makes the optimal adjustment of the controller more difficult.

[0011] When a PI or PID controller is used, the degrees of the associated polynomials are two and three coefficients, respectively. The adjustment is therefore simple. However, to achieve the objectives of robustness in performance and stability is no longer possible in the presence of an electromechanical drive train characterized by resonant modes, and possibly by non-linearities such as transmission play.

[0012] To remedy the difficulties that occur when a PI or PID controller is used in an electromechanical drive train with little damping, it has been proposed to associate with said controller a state feedback based on a signal issued by an estimator. This solution is described in the thesis of Marius Goslar entitled “Ein Beitrag zur anwendungsorientierten Zustands-regelung elektrischer Hochleistungsantriebe,” presented on Aug. 14, 1998 and published in Great Britain by “Conservatree Print & Design” ISBN 0953473503.

[0013] More precisely, the abovementioned thesis concerns a case in which the angular speed of the load is known. It is proposed to use the genetic algorithms to optimize the coefficients of the polynomials of the controller.

[0014] This known solution provides good damping of the system when resonance frequencies are present. However, it is still not very satisfactory in terms of robustness. Moreover, the use of the genetic algorithms to calculate the coefficients of the polynomials of the controller makes the adjustment of said controller particularly complex.

[0015] When an H_(∞), QFT, or RST controller is used, the degrees of the polynomials are appreciably higher. These controllers are called “robust” because they allow the controlled system to preserve its stability and performance under nominal conditions, that is, where there are uncertainties about the parameters, and possibly nonlinearities in the transmission, such as transmission play. However, they have a high number of coefficients that makes their adjustment particularly complex.

[0016] Regardless of the type of controller used, the robust control of the angular speed of the load is made difficult when it is not possible to have a direct measurement of this speed.

[0017] A specific object of the invention is a method and device for controlling the mechanical speed of a load driven in rotation through an electromechanical drive train with little damping, without measuring this speed, while presenting satisfactory nominal static and dynamic performances and good robustness, for a minimum number of control coefficients.

[0018] According to the invention, this result is achieved by a method of controlling the speed of a load driven in rotation by an electric motor through an electromechanical drive train with little damping, according to which the real speed of the motor is measured and a torque actuation reference quantity Γ* of the motor with the measured speed Ω_(m) of the motor and a reference speed Ω_(L)* of the load are determined, characterized in that said reference quantity Γ* is determined:

[0019] by calculating an estimated speed {circumflex over (Ω)}_(L) of the load by an observer the inputs of which are the reference quantity Γ* and an innovation error represented by the weighted difference between the measured speed Ω_(m) of the motor and an estimated speed Ω_(m) of the motor;

[0020] by calculating the difference between the reference speed Ω_(L)* of the load and the estimated speed {circumflex over (Ω)}_(L) of the load; and

[0021] by using this calculated difference as input of a robust CRONE-type non-integer order controller, the physical realization of which is provided by an integer order controller, obtained by approximation in the frequency domain.

[0022] In the method thus defined, the calculation of the estimated speed of the load enables an effective control of this speed in the absence of its measurement. In practice, this calculation is done by an observer, such as a Kalman filter, suitable for reconstructing the internal state of the electromechanical drive train from a dynamic model thereof.

[0023] Appropriately, the observer uses a dynamic model such as: ${\begin{bmatrix} {\overset{.}{X}}_{CC} \\ {\overset{.}{X}}_{AE} \\ {\overset{.}{\Gamma}}_{L} \end{bmatrix} = {{\begin{bmatrix} \left\lbrack A_{CC} \right\rbrack & \begin{bmatrix} {couplage} \\ {1/2} \end{bmatrix} & \begin{bmatrix} {couplage} \\ {{let2}/3} \end{bmatrix} \\ \lbrack 0\rbrack & \left\lbrack A_{AE} \right\rbrack & \begin{bmatrix} {couplage} \\ {{let2}/3} \end{bmatrix} \\ \lbrack 0\rbrack & \lbrack 0\rbrack & \lbrack 0\rbrack \end{bmatrix}\begin{bmatrix} X_{CC} \\ X_{AE} \\ \Gamma_{L} \end{bmatrix}} + {\begin{bmatrix} B \\ 0 \end{bmatrix}\Gamma^{*}}}},\left( {{couplage} = {coupling}} \right)$

[0024] in which A_(CC), A_(AE), and B are matrices of coefficients assigned respectively to the drive train, to the electrical actuator, and to the nature of the disturbance of the resisting torque applied to the load, X_(CC), X_(AE), and Γ_(L) representing respectively the status of the drive train, the electrical actuator and the resisting torque applied to the load. Moreover, the expressions “couplage1/2” and “couplagelet2/3” designate respectively a matrix of coupling coefficients between the electrical actuator and the drive train and a matrix of coupling coefficients between the electrical actuator-drive train unit and the load.

[0025] Furthermore, the use of a non-integer order controller to produce the integer order controller being used to correct the estimated speed of the load makes it possible to benefit from the properties of robustness of this type of control while limiting the number of control coefficients to a value comparable to that of non-robust PI or PID type controls. Thus a robust performance is ensured of the speed control with respect to variations of the parameters of the drive train, such as inertia and rigidity, and with respect to the non-linearities of the system such as transmission play. In practice, the number of control coefficients here is equal to three, which considerably simplifies the adjustment of the controller.

[0026] These advantages derive from the remarkable properties of a robust non-integer order control, or CRONE (registered trademark). The characteristics of this control are described in detail in the work “La Commande CRONE, du Scalaire au Multivariable” [CRONE control, from Scalar to Multivariable] by A. Oustaloup and B. Mathieu, Hermes, 1999.

[0027] Appropriately, and according to the characteristics peculiar to the third generation CRONE control as set forth in the above-mentioned work, a non-integer order control is used, the transmittance β(s) of which is such that: ${{\beta (s)} = {y_{0}\left\lbrack \left( \frac{\omega_{s}}{s} \right)^{n} \right\rbrack}_{C_{j}}},{{{with}\quad n} \in C_{j}},$

[0028] where n=a+ib is a complex non-integer order of integration of which the real part a determines the phase placement of the pattern at the central frequency ω_(g) and of which the imaginary part b determines the slope of the pattern with reference to the vertical, C_(j) in the Nichols plan, and y₀ is a gain placement. The letters a, b, and ω_(g) correspond to the controller's adjustment coefficients. The central frequency ω_(g) is relative to the controller's response time. It is also used in the relationships $\omega_{A} = \frac{\omega_{g}}{k}$

[0029] and ω_(B)=kω_(g), in which ω_(A) and ω_(B) are extreme frequencies defining the range of frequencies in which the robustness of performance and stability should be ensured, and k is a dependent coefficient of said range of frequencies.

[0030] In the case of an electromechanical drive train having multiple resonances, as many rejection filters as there are resonant modes are associated with this transmittance, as described in section 3.4 (“Cas des procédés résonants” [cases of resonant processes]) in the abovementioned work.

[0031] Preferably, the integer order controller is obtained by effecting a frequency truncation at low and high frequencies, defined by low ω_(b) and high ω_(h) frequency boundaries, and by calculating, by approximation in this range of frequencies (ω_(b) ω_(h)), the equivalent integer order transmittance C_(N)(S) from the following formula: ${{C_{N}(s)} = {C_{0}{\prod\limits_{k = {- N}}^{N}\frac{1 + \frac{s}{{\omega^{\prime}}_{k}}}{1 + \frac{s}{\omega_{k}}}}}},$

[0032] in which: ${\frac{{\omega^{\prime}}_{k + 1}}{{\omega^{\prime}}_{k}} = {\frac{\omega_{k + 1}}{\omega_{k}} + {\alpha \quad \eta}}},{\frac{{\omega^{\prime}}_{k + 1}}{\omega_{k}} = {{\eta \quad {and}\quad \frac{\omega_{k}}{{\omega^{\prime}}_{k}}} = \alpha}},$

[0033] with α>1 and η>1, and where C₀ is the gain at zero frequency.

[0034] The invention also concerns a device for controlling the speed of a load driven in rotation by an electric motor through an electromechanical drive train with little damping, comprising means for measuring the real speed of the motor and means for determining a torque actuation reference quantity Γ* from a speed signal Ω_(m) issued by said speed measuring means and from a set point signal Ω_(L)* representative of a reference speed of the load, characterized in that the means for determining said reference quantity include:

[0035] an observer the inputs of which are the reference quantity Γ* and an innovation error represented by the weighted difference between the measured speed Ω_(m) of the motor and an estimated speed {circumflex over (Ω)}_(m) of the motor, the observer being able to calculate an estimated speed {circumflex over (Ω)}_(L) of the load and the difference between the reference speed Ω_(L)* of the load and the estimated speed {circumflex over (Ω)}_(L) of the load;

[0036] a robust CRONE-type non-integer order controller, using as an input the calculated difference between the reference speed Ω_(L)* of the load and the estimated speed {circumflex over (Ω)}_(L) of the load and the physical realization of which is provided by a integer order controller, obtained by approximation in the frequency domain.

[0037] By way of non-limiting example, a preferred form of embodiment of the invention will now be described, with reference to the attached drawings in which:

[0038]FIG. 1, already described, diagrammatically represents a closed loop servo system according to the prior art;

[0039]FIG. 2 is a view comparable to FIG. 1, which diagrammatically represents a servo system that includes a device for controlling the angular speed of the load according to the invention;

[0040]FIG. 3 is a block diagram that represents in more detail the CRONE controller of the control device according to the invention;

[0041]FIG. 4 represents very diagrammatically a model of an electromechanical drive train used to test the control device according to the invention;

[0042]FIG. 5 represents, at A, the development curve of the speed Ω_(L) of the load driven by the electromechanical drive train of FIG. 4 (in reduced units) as a function of time t (in seconds) and, at B, the development curve of the actuating torque Γ* (in reduced units) delivered by the controller as a function of time t, when a PI controller of the prior art and a system like the one illustrated in FIG. 1 are used;

[0043]FIG. 6 represents, at A′, the development curve of the speed Ω_(L) of the load driven by the electromechanical drive train of FIG. 4 (in reduced units) issued by the controller as a function of time t (in seconds) and, at B′, the development curve of the actuating torque Γ* (in reduced units) delivered by the controller as a function of time t, when an observer (Kalman filter) is associated with a non-integer CRONE controller, according to the invention and according to the assembly illustrated in FIG. 2;

[0044]FIG. 7 represents the transfer functions (amplitude of the signal in dB, as a function of the frequency in Hz) of the drive train corresponding to the model of FIG. 4, respectively when the inertias have their nominal values (curve C) and in two extreme cases of variation of said inertias (curves D and E); and

[0045]FIG. 8 represents the curves A′ and B′ of FIG. 6 in the case when the inertias have their nominal values and in two extreme cases of variation of said inertias, on the model of FIG. 4.

[0046] As illustrated very diagrammatically in FIG. 2, a speed control device according to the invention is designed to control the speed of rotation of a load 10 driven in rotation by means of an electromechanical drive train. Said electromechanical drive train includes a torque actuator 12 that drives the load 10 in rotation through a drive train 13. The torque actuator 12 has an electric motor 14 coupled to a power converter 16.

[0047] It is assumed that the electromechanical drive train has little damping and there are no means to measure the mechanical speed of the load 10.

[0048] According to the invention, the control of the speed of rotation of the load 10 is provided by associating an observer 18 and a CRONE-type non-integer order controller 20.

[0049] More precisely, a function of the observer 18 is to calculate an estimated speed of rotation {circumflex over (Ω)}_(L) of the load 10. To that end, the observer 18 receives at a first input a first signal Γ*, representative of a reference quantity of torque actuation of the motor 14. The observer 18 calculates at 22 an innovation error, represented by the difference between the real angular speed Ω_(m) of the motor 14 and an estimated angular speed {circumflex over (Ω)}_(m) of the motor 14. This difference is weighted by a weighting coefficient K, as represented at 24 in FIG. 2.

[0050] The reference quantity of the motor's actuation torque, represented by the first signal Γ*, corresponds to a desired actuation torque of the motor 14. The first signal Γ* is taken directly at the output of the CRONE controller 20 before being injected into the torque actuator 12.

[0051] For purposes of simplification, the speed of rotation of the output shaft of the motor 14 is called “real angular speed Ω_(m) of the motor.” This speed can be measured by any known means such as a speed sensor 26, without going beyond the scope of the invention. A signal representing the real angular speed Ω_(m) of the motor is received at a second input of the observer 18.

[0052] The observer 18 incorporates a dynamic model 28 that calculates the estimated speed of rotation {circumflex over (Ω)}_(L) of the load 10 as well as the estimated angular speed {circumflex over (Ω)}_(m) of the motor 14, used to calculate the innovation error. To that end, the observer 18 reconstructs the electromechanical drive train and continuously resets the states of the system by calculating the innovation error at 22.

[0053] This dynamic model rests on the internal representation of the system, usually expressed by the equations: $\left\{ {\begin{matrix} {\overset{.}{\hat{X}} = {{A\hat{X}} + {BU} + {K\left( {\hat{Y} - Y} \right)}}} \\ {\hat{Y} = {{C\hat{X}} + {DU}}} \end{matrix},} \right.$

[0054] In these equations, U represents the vector of the inputs from the electromechanical drive train, X represents an element of the space state of the system, Y represents the vector of the outputs of the system, and A, B, C, and D are matrices the coefficients of which are real numbers. The symbols X and Y represent the estimated values of X and Y obtained by the implementation of the model.

[0055] In practice, the observer 18 can be comprised of a Kalman filter. In this case, the calculations are done linearly and in discrete time. In particular, this type of observer has the advantage of having minimal variations. However, any other observer that can obtain comparable results, such as a Luenberger filter, can be used without going beyond the scope of the invention.

[0056] In order to ensure a faithful representation of the internal state of the system, the model takes into account the electromagnetic model and the mechanical model of the electromechanical drive train, as well as the disturbances occurring on the resisting torque applied to the load.

[0057] Thus, the representation of the torque actuator 12 comprising the power converter 16 and the motor 14 is provided by a first continuous equivalent state model, the free development of which is represented by:

[{dot over (X)} _(AE) ]=[A _(AE) ][X _(AE)].

[0058] In a comparable way, the drive train 13 is represented by a second continuous equivalent state model, the free development of which is represented by:

8 {dot over (X)} _(CC) ]=[A _(CC) ][X _(CC)].

[0059] The interconnection of the two units thus defined leads to representing the electromechanical drive train in a form the free development of which is represented by: $\begin{bmatrix} {\overset{.}{X}}_{CC} \\ {\overset{.}{X}}_{AE} \end{bmatrix} = {{\begin{bmatrix} \left\lbrack A_{CC} \right\rbrack & \begin{bmatrix} {couplage} \\ {1/2} \end{bmatrix} \\ \lbrack 0\rbrack & \left\lbrack A_{AE} \right\rbrack \end{bmatrix}\begin{bmatrix} X_{CC} \\ X_{AE} \end{bmatrix}}.}$

[0060] In this relationship, the expression “couplage1/2” [coupling1/2] designates a matrix of coupling coefficients between the torque actuator 12 and the drive train 13.

[0061] In order also to take into account the disturbances of the torque Γ_(L) applied to the load, it is assumed that {dot over (Γ)}_(L)=0.

[0062] The system can then be stated in the form of a dynamic model so that the response to the torque setting input Γ* can be represented by: $\begin{bmatrix} {\overset{.}{X}}_{CC} \\ {\overset{.}{X}}_{AE} \\ {\overset{.}{\Gamma}}_{L} \end{bmatrix} = {{\begin{bmatrix} \left\lbrack A_{CC} \right\rbrack & \begin{bmatrix} {couplage} \\ {1/2} \end{bmatrix} & \begin{bmatrix} {couplage} \\ {{let2}/3} \end{bmatrix} \\ \lbrack 0\rbrack & \left\lbrack A_{AE} \right\rbrack & \begin{bmatrix} {couplage} \\ {{let2}/3} \end{bmatrix} \\ \lbrack 0\rbrack & \lbrack 0\rbrack & \lbrack 0\rbrack \end{bmatrix}\begin{bmatrix} X_{CC} \\ X_{AE} \\ \Gamma_{L} \end{bmatrix}} + {\begin{bmatrix} B \\ 0 \end{bmatrix}{\Gamma^{*}.}}}$

[0063] In this relationship, the expression “couplagelet2/3” [coupling1 and2/3] designates a matrix of coupling coefficients between the unit formed by the torque actuator 12 and the drive train 13, and the load 10. The coefficients of the coupling matrices “couplage1/2” and “couplagelet2/3” as well as the different coefficients of the matrices A_(CC), A_(AE), and B are defined by a person skilled in the art as a function of the respective characteristics of the drive train 13, the torque actuator 12, and the disturbances of the resisting torque applied to the load 10. As can be seen, this operation derives from the simple implementation of basic knowledge of a specialist in automation.

[0064] As is diagrammatically represented in FIG. 3, the controller 20 is comprised of a robust CRONE-type non-integer order controller 32, the physical realization of which is provided by an integer order controller 34, obtained by approximation in the frequency domain.

[0065] The robust non-integer order controller 32 is a third generation CRONE-type controller, that is, it uses a generalized non-integer order frequency pattern. The transmittance β(s) of this pattern is given by the relationship: $\begin{matrix} {{{\beta (s)} = {y_{0}\left\lbrack \left( \frac{\omega_{s}}{s} \right)^{n} \right\rbrack}_{C_{j}}},{{{with}\quad n} \in C_{j}},} & (1) \end{matrix}$

[0066] in which n is a complex non-integer order of integration, such as n=a+ib. The coefficients a and b as well as the term ω_(g) constitute the only adjustment coefficients of the frequency pattern used by the robust controller 32. In other words, the adjustment of the controller is made as soon as the values of these three coefficients are established. This characteristic constitutes an essential advantage of the invention, because it enables an adjustment that is practically as simple as in the case when a non-robust PI or PID controller is used.

[0067] The real part a of n determines the placement in phase of the pattern at the central frequency ω_(g), in the Black or Nichols plan, designated by C_(j). The imaginary part b determines the slope of the pattern with reference to the vertical in said plan.

[0068] In the relationship (1) above, y₀ represents a gain placement. The central frequency ω_(g) is relative to the controller's response time. It is also used to determine the frequencies $\omega_{A} = \frac{\omega_{g}}{k}$

[0069] and ω_(B)=kω_(g), with ωε[ω_(A),ω_(B)]. In other words, the letters ω_(A) and ω_(B) designate the extreme frequencies defining the range of frequencies in which the robustness in performance and stability can be assured. Finally, the letter k designates a coefficient that depends on said range of frequencies defined between ω_(A) and ω_(B).

[0070] For more details concerning the robust controller 32 and a summary thereof, it will be useful to refer to the abovementioned work of Alain Oustaloup and Benôit Mathieu, and particularly Chapter 2 entitled: “Commande CRONE de troisième génération” [third generation CRONE control].

[0071] When the electromechanical drive train has multiple resonances, as many rejection filters are associated with the transmittance β(s) of the frequency pattern as there are resonant modes in said drive train. It will be useful to refer to section 3.4 of the abovementioned work of Alain Oustaloup and Benôit Mathieu, incorporated here by reference.

[0072] The integer order controller 34 comprises a robust monovariable controller. It functions in discrete time with a sampling period that can be greater than or equal to that of the observer 18. Thus, in the form of embodiment described here by way of example, the sampling period of the integer order controller 34 is five times longer than that of the observer 18.

[0073] The integer order controller 34 receives a representative signal of the difference (number 30 in FIG. 2) between the reference speed Ω_(L)* of the load 10 and the estimated speed {circumflex over (Ω)}_(L) of the load. The reference speed Ω_(L)* of the load 10 represents the desired or set speed of rotation of the load. The estimated speed {circumflex over (Ω)}_(L) of the load 10 is furnished by the observer 18. Based on this difference, the integer order controller issues a representative signal Γ* of the torque that should be applied to the motor 14 by the power converter 16.

[0074] The determination of the discrete integer order controller 34 is done by first performing, if necessary, a low and high frequency truncation of the non-integer order controller 32, by means of an integer order transmittance β_(b)(s,ω_(b)) defined by a high frequency boundary ω_(b) and an integer order transmittance β_(h)(s,ω_(h)) defined by a high frequency boundary ω_(h). These two frequency boundaries can be symmetrically distributed with reference to the extreme frequencies ω_(A) and ω_(B) of the generalized pattern, so that $\omega_{b} = \frac{\omega_{A}}{K^{\prime}}$

[0075] and ω_(h)=K′ω_(B), K′ being a coefficient equal to 10, for example, so that $\frac{\omega_{g}}{\omega_{b}}\operatorname{>>}{1\quad {and}\quad \frac{\omega_{g}}{\omega_{h}}{\operatorname{<<}1.}}$

[0076] This frequency truncation makes it possible to desensitize the controller 32 with respect to the high frequencies (e.g., measurement noises) and to ensure a zero static error. For more details on this subject, it will be useful to refer to the work of Alain Oustaloup and Benoit Mathieu as cited above (see in particular the Bode diagram of FIG. 1.3). The calculation of the adjustment coefficients of the controller is then done on the basis of the frequency truncated non-integer order controller 32′ (FIG. 3).

[0077] The integer order controller 36 is then determined by an approximation in the frequency domain of the non-integer controller 32′. The transmittance C_(N)(S) of the integer order controller 36 is then given by the relationship: $\begin{matrix} {{{C_{N}(s)} = {C_{0}{\prod\limits_{k = {- N}}^{N}\frac{1 + \frac{s}{{\omega^{\prime}}_{k}}}{1 + \frac{s}{\omega_{k}}}}}},} & (2) \end{matrix}$

[0078] in which: ${\frac{{\omega^{\prime}}_{k + 1}}{{\omega^{\prime}}_{k}} = {\frac{\omega_{k + 1}}{\omega_{k}} = {\alpha\eta}}},{\frac{{\omega^{\prime}}_{k + 1}}{\omega_{k}} = {{\eta \quad {et}\quad \frac{\omega_{k}}{{\omega^{\prime}}_{k}}} = \alpha}},$

[0079] with α>1 and η>1 and where C₀ represents the zero frequency gain.

[0080] The integer order controller 36 is then discretized on the basis of a sampling frequency that can be greater than or equal to that of the observer 18, in order to obtain the discrete integer order controller 34 that will be used in the real time calculator.

[0081] The observer 18 as well as the controllers 32, 32′, 34, and 36 use one or more calculators.

[0082] Represented in FIG. 4 is a model that was made to validate through testing, in a concrete case, the results obtained by the implementation of the speed control method according to the invention. More specifically, comparative tests were made on the model of FIG. 4, using successively a PI controller of the prior art and a CRONE type observer controller unit according to the invention.

[0083] The model used for the tests comprised a first motor 114, “speed control,” comparable to the motor 14 in FIG. 2. The motor 114 was mechanically connected to a second motor 110, “the load,” through a drive train 113 comprised of a first extension piece 113A, a first inertial test weight 113B, a second extension piece 113B, a second inertial test weight 113D, a third extension piece 113E that includes a spacer 113F and a third inertial test weight 113G.

[0084] In the model thus produced, the motor 110 acts as the load 10 in FIG. 2 and the drive train 113 acts as the train 13.

[0085] Each of the motors 114 and 110 had an inertia of 0.00488 m²*kg and the rigidity of the extension pieces 113A, 113C, and 113E were respectively 357 N*m/rad, 175 N*m/rad, and 78 N*m/rad. Finally, the nominal inertias of the inertial weights 113B, 113D, and 113G were respectively equal to 0.014 m²*kg, 0.03 m²*kg, and 0.068 m²*kg, while still being able to simulate large errors causing these inertias to vary between one half and twice their nominal values.

[0086] The tests performed on this model by using a PI controller of the prior art are illustrated by the curves A and B in FIG. 5. More specifically, the curve A shows the development of the speed of rotation Ω_(L) of the load (in reduced units) as a function of the time t (in seconds), and the curve B shows the development of the torque Γ* delivered by the controller (in reduced units) as a function of the time t (in seconds). As can be seen, a disturbance of the load torque at the moment t₀ results in particular in a relatively long period before return to equilibrium and by the system entering into resonance.

[0087] As curves A′ and B′ in FIG. 6 illustrate, the situation is very different when the speed control device according to the invention is used. More precisely, the curve A′ shows the development of the speed of rotation Ω_(L) of the load (in reduced units) as a function of the time t (in seconds) and the curve B′ shows the development of the torque Γ* delivered by the controller (in reduced units) as a function of the time t (in seconds). As can be seen, in this case the system very quickly returns to equilibrium after the disturbance of the load torque at the time t₀. Moreover, no resonance phenomenon occurs, unlike in the prior art.

[0088] The curves in FIG. 6 show that the speed control device according to the invention has the desired damping capacities when there are resonance frequencies in the electromechanical drive train.

[0089] Another series of tests was performed to show the robustness capacity of the speed control device according to the invention. These tests are illustrated in FIGS. 7 and 8.

[0090] Represented in FIG. 7 is the transfer function of the electromechanical drive train of the model illustrated in FIG. 4, that is, the development of the amplitude of the signal (in dB) as a function of the frequency (in Hz), respectively in the case where the inertial weights 113B, 113D, and 113G are nominal (curve C), in the case where the inertias of the first weight 113B and third weight 113G are divided by two and second weight 113D is multiplied by two (curve D), and in the case where the intertias of the first weight 113B and third weight 113G are divided by two and the second weight 113D is nominal (curve E). As can be seen, the resonance frequencies of the electromechanical drive train are significantly modified.

[0091] In FIG. 8, the curves A′ and B′ are shown that were obtained in the cases corresponding to the transfer functions illustrated at C, D, and E of FIG. 7. It can be seen that, in spite of the large variations in the inertias and in the resonance frequencies resulting therefore, the closed loop system remains stable.

[0092] These experiments confirm that the speed control according to the invention has excellent properties of robustness in spite of the absence of measurement of the controlled speed, unlike the methods of the prior art using PI controllers.

[0093] Furthermore, the association of an observer and a CRONE-type controller, according to the invention, makes it possible to obtain these properties of robustness with only three adjustment coefficients (the real a and imaginary b parts of the complex non-integer integration order n and the frequency ω_(g)).

[0094] Of course, the method and device according to the invention can be applied to any system comprising a load driven in rotation by an electromechanical drive train with little damping, regardless of the nature of the load and irrespective of the characteristics of the electromechanical drive train, and particularly its resonance frequencies and possible play. 

1. Method of controlling the speed of a load (10) driven in rotation by an electric motor (12) through an electromechanical drive train (13) with little damping, according to which the real speed of the motor (12) is measured and a torque actuation reference quantity Γ* of the motor (12) with a measurement speed Ω_(m) of the motor and a reference speed Ω_(L)* of the load (10) are determined, characterized in that said reference quantity Γ* is determined: by calculating an estimated speed {circumflex over (Ω)}_(L) of the load by an observer (18) the inputs of which are the reference quantity Γ* and an innovation error represented by the weighted difference between the measured speed Ω_(m) of the motor and an estimated speed {circumflex over (Ω)}_(m) of the motor; by calculating the difference between the reference speed Ω_(L) *of the load (10) and the estimated speed {circumflex over (Ω)}_(L) of the load; and by using this calculated difference as input of a robust CRONE-type non-integer order controller (20), the physical realization of which is provided by an integer order controller, obtained by approximation in the frequency domain.
 2. Method according to claim 1, in which the estimated speed {circumflex over (Ω)}_(L) of the load (10) is calculated from a dynamic model of the electromechanical drive train (13).
 3. Method according to claim 2, in which a dynamic model is used that takes into account the electromagnetic model, the mechanical model, and the disturbances of the torque of the load (10).
 4. Method according to claim 3, in which is used a dynamic model such as: $\begin{matrix} {\begin{bmatrix} {\overset{.}{X}}_{CC} \\ {\overset{.}{X}}_{AE} \\ {\overset{.}{\Gamma}}_{L} \end{bmatrix} = \begin{bmatrix} \left\lbrack A_{CC} \right\rbrack & \left\lbrack {{couplage1}/2} \right\rbrack & \left\lbrack {{couplage1et2}/3} \right\rbrack \\ \lbrack 0\rbrack & \left\lbrack A_{AE} \right\rbrack & \left\lbrack {{couplage1et2}/3} \right\rbrack \\ \lbrack 0\rbrack & \lbrack 0\rbrack & \lbrack 0\rbrack \end{bmatrix}} \\ {{{\begin{bmatrix} X_{CC} \\ X_{AE} \\ \Gamma_{L} \end{bmatrix} + {\begin{bmatrix} B \\ 0 \end{bmatrix}\Gamma^{*}}},}} \end{matrix}$

which A_(CC), A_(AE), and B are matrices of coefficients assigned respectively to the drive train (16), to the electrical actuator (12, 14), and to the nature of the disturbance of the resisting torque applied to the load (10), X_(CC), X_(AE), and Γ_(L) representing respectively the status of the drive train, the electrical actuator, and the resisting torque applied to the load. Moreover, the expressions “couplage1/2” and “couplagelet2/3,” designate respectively a matrix of coupling coefficients between the electrical actuator and the drive train and a matrix of coupling coefficients between the electrical actuator-drive train unit and the load.
 5. Method according to any of claims 2 to 4, in which a Kalman filter is used in discrete time as observer (18).
 6. Method according to any of claims 1 to 5, in which a non-integer order controller is used the transmittance β(s) of which is such that: ${{\beta (s)} = {y_{0}\left\lbrack \left( \frac{\omega_{g}}{s} \right)^{n} \right\rbrack}_{C_{j}}},$

with nεC_(j), where n=a+ib is a complex non-integer order of integration of which the real part a determines the phase placement of the pattern at the central frequency ω_(g) and of which the imaginary part b determines the slope of the pattern with reference to the vertical, C_(j) in the Nichols plan, and y₀ is a gain placement and the frequency ω_(g) is relative to the controller's response time and is also used in the relationships $\omega_{A} = \frac{\omega_{g}}{k}$

and ω_(B)=kω_(g), in which ω_(A) and ω_(B) are extreme frequencies defining the range of frequencies in which the robustness of performance and stability should be ensured, and k is a dependent coefficient of said range of frequencies.
 7. Method according to claim 6, in which the electromechanical drive train (13) having multiple resonances, as many rejection filters are associated with the transmittance β(s) of the frequency pattern as there are resonant modes in said drive train.
 8. Method according to either of claims 6 and 7, in which the integer order controller is obtained by effecting a frequency truncation at low and high frequencies, defined by low ω_(b) and high ω_(h) frequency boundaries, and by calculating, by approximation in this range of frequencies (ω_(b) ω_(h)), the equivalent integer order transmittance C_(N)(s) from the following formula: ${{C_{N}(s)} = {C_{0}{\prod\limits_{k = {- N}}^{N}\frac{1 + \frac{s}{{\omega^{\prime}}_{k}}}{1 + \frac{s}{\omega_{k}}}}}},$

in which: ${\frac{{\omega^{\prime}}_{k + 1}}{{\omega^{\prime}}_{k}} = {\frac{\omega_{k + 1}}{\omega_{k}} = {\alpha\eta}}},{\frac{{\omega^{\prime}}_{k + 1}}{\omega_{k}} = {{\eta \quad {and}\quad \frac{\omega_{k}}{{\omega^{\prime}}_{k}}} = \alpha}},$

with α>1 and η>1, and where C₀ is the gain at zero frequency.
 9. Device for controlling the speed of a load (10) driven in rotation by an electric motor (12) through an electromechanical drive train (13) with little damping, comprising means (26) for measuring the real speed of the motor (12) and means for determining a torque actuation reference quantity Γ* of the motor (12) from a speed signal Ω_(m) issued by said speed measuring means (26) and from a set point signal Ω_(L)* representative of a reference speed of the load (10), characterized in that the means for determining said reference quantity include: an observer (18) the inputs of which are the reference quantity Γ* and an innovation error represented by the weighted difference between the measured speed Ω_(m) of the motor (12) and an estimated speed {circumflex over (Ω)}_(m) of the motor, the observer (18) being able to calculate an estimated speed {circumflex over (Ω)}_(L) of the load and the difference between the reference speed Ω_(L)* of the load (10) and the estimated speed Ω_(L) of the load; a robust CRONE-type non-integer order controller (20), using as an input the calculated difference between the reference speed Ω_(L)* of the load (10) and the estimated speed {circumflex over (Ω)}_(L) of the load and the physical realization of which is provided by a integer order controller, obtained by approximation in the frequency domain.
 10. Device according to claim 9, in which the observer (18) calculates the estimated speed {circumflex over (Ω)}_(L) of the load (10) from a dynamic model of the electromechanical drive train (13).
 11. Device according to claim 10, in which the dynamic model includes three interconnected matrix units, taking into account the electromagnetic model, the mechanical model, and the disturbances of the torque of the load (10), respectively.
 12. Device according to any of claims 9 to 11, in which the observer (18) is a Kalman filter in discrete time.
 13. Device according to any of claims 9 to 12, in which the non-integer order controller (20) has a non-integer order transmittance β(s) such that: ${{\beta (s)} = {y_{0}\left\lbrack \left( \frac{\omega_{g}}{s} \right)^{n} \right\rbrack}_{C_{j}}},$

with nεC_(j), where n=a+ib is a complex non-integer order of integration of which the real part a determines the phase placement of the pattern at the central frequency ω_(g) and of which the imaginary part b determines the slope of the pattern with reference to the vertical, C_(j) in the Nichols plan, and y₀ is a gain placement and the frequency ω_(g) is relative to the controller's response time and is also used in the relationships $\omega_{A} = \frac{\omega_{g}}{k}$

and ω_(B)=kω_(g), in which ω_(A) and ω_(B) are extreme frequencies defining the range of frequencies in which the robustness of performance and stability should be ensured, and k is a dependent coefficient of said range of frequencies.
 14. Device according to claim 13, in which the electromechanical drive train has multiple resonances and the controller (20) has as many rejection filters as there are resonant modes in said drive train.
 15. Device according to either of claims 13 and 14, in which the integer order controller effects a frequency truncation at low and high frequencies, defined by low ω_(b) and high ω_(h) frequency boundaries, and has an integer order transmittance C_(N)(S) such that: ${{C_{N}(s)} = {C_{0}{\prod\limits_{k = {- N}}^{N}\quad \frac{1 + \frac{s}{\omega_{k}^{\prime}}}{1 + \frac{s}{\omega_{k}}}}}},$

in which: $\begin{matrix} {{\frac{\omega_{k + 1}^{\prime}}{\omega_{k}^{\prime}} = {\frac{\omega_{k + 1}}{\omega_{k}} = {\alpha \quad \eta}}},} \\ {{\frac{\omega_{k + 1}^{\prime}}{\omega_{k}} = {{\eta \quad {and}\quad \frac{\omega_{k}}{\omega_{k}^{\prime}}} = \alpha}},} \end{matrix}$

with α>1 and η>1, and where C₀ is the gain at zero frequency. 